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We consider an uniformly charged incompressible nuclear fluid bounded by a closed surface. It is 
shown that an evolution of an axisymmetric surface T(r, t) = a — E(z, t) = 0, r — (a, <f), z) can be 
approximately reduced to a motion of a curve in the (a, z)-plane. A nonlinear integro-diffrerential 
equation for the contour T,(z,t) is derived. It is pointed on a direct correspondence between £(jz,t) 
and a local curvature, that gives possibility to use methods of differential geometry to analyze an 
' evolution of an axisymmetric nuclear surface. 
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I. MOTIVATION 



Nonlinear dynamics of a nuclear surface is an object of special interest due to the following reasons. First, nuclear 
density falls considerably in the surface region, where the density fluctuations and clustering may be important. 
Second, different types of instability may be developed in the surface region and lead to fragmentation processes 
at low (fission, nucleon transfer) and high (multifragmentation, break-up etc.) energies. Finally, the liquid drop 
£S| ' concept |j| for over a century are intensively used in macro- Q and micro-physics || . Nonlinear dynamics of shapes 
7-H , in any complicated systems enevitably leads to mathematical problem of describing global geometric quantities such 
1 as surface and enclosed volume in different dimensions (polymers, cell membranes, 3D droplets). 

An application of a soliton concept to nonlinear nuclear hydrodynamics has given yet new possibilities in this field 
(See e.g. review [Q and for the recent refs.). However any extension of nonlinear dynamics from 1+1 to 2+1 and 

3+1 dimensions meets principal difficulties. The crucial point is to reduce the dimension of a problem. In paper j^] 
we considered the simplest two-dimensional nonlinear liquid objects. It was shown that 2D pure vortical motion of 
inviscid nuclear liquid can be reduced to ID evolution of the contour bounding this drop. In the next paper fie] ] the 
extension to semi-3D geometry was done. The equations of motion describing localized vortex on a spherical nuclear 
surface - a bounded region of constant vorticity, surrounded by irrotational flux - were reduced to ID nonlinear 
evolution of the boundary. 

In this short report we consider an uniformly charged incompressible nuclear 3D fluid bounded by a closed surface. 
It is shown that evolution of an axisymmetric surface T(r, t) = a — £(z, i) = 0, r = (c, <f>, z) can be approximately 
reduced to a motion of a curve in the (cr, z)-plane. 



II. FRAMEWORK 



The evolution of one-body Wigner phase-space distribution function is analyzed instead of a full many-body wave 
function. Integrating the kinetic equation 

df p df dV df _ p 2 

dt m dr dr dp rei ' W 2m 

over the momentum space with different polynomial weighting functions of the p-variable one comes to an infinite 
chain of equations for local collective observables including the density, collective velocity, pressure and an infinite set 
of tensorial functions of the time and space coordinates, which are defined as moments of the distribution function in 
the momentum space: 

- the particle n(r,t) = g J dp f(v,p,t), and the mass p{r,t) = mn(r,t) densities, where we consider a proton and 
neutron as different state of identical particle (the spin-isospin degeneration g = 4) . 

- the collective current and velocity of nuclear matter p(r,t)u(r,t) = g J dppf(r,p,t) 
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the pressure tensor Py(r,i) = [g/m) f dpqiqjf(r,p,t), q = p — mu 
the energy and momentum transfer tensors of a different order 

Pijk(r,t) = J d Pqiqjqkf(r,p,t), P^ k (r,t) = J dpqiqj..qkf(r,p,i), (2) 

and the integrals related to relaxation terms 

dpl re i = 0, J dpplrei = 0, Rij = £- J dpqiqjl reh ... (3) 
Truncating this chain one arrives at the "fluid dynamical" level of description of nuclear processes. 



Dili . s^dP lk p dV /n x^ n (V i , 

k k 

+ p y ^2 / £ isi (2Q, s Uj + — -Xj) = 0, (5) 
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where the usual notation = J| + X) u fc af^ ^ s introduced for the operator giving the material derivative, or rate 
of change at a point moving with the fluid locally. The hydrodynamical set of Eqs. (||H) describes an evolution in a 

3 

rotating nuclear system. The linear transformation Xi = X TijXj relates the coordinates of a point, (X\, X2, X3), in 

3=1 

an inertial frame, and,(xi, X2, £3 in a moving frame of reference with a common origin. The orientation of a moving 
frame, with respect to an inertial frame, will be assumed to be time dependent, Ty(i), representing an orthogonal 
transformation, and the vector 

°i = 5 E E »k{dT/dt). m T+ k 

represents a general time-dependent rotation. 

Let us restrict ourselves with simplest possible motion of nuclear fluid. Consider an uniformly charged (with a 
total charge Ze) incompressible nuclear "fluid" confined to a volume V bounded by a closed surface S, obeying 
to the equation T(r,t) = a — T,(z,t) = 0, where (er, 4>, z) are cylindrical coordinates of a point. The geometry 
places the origin at the center-of-mass (/ dr p(r , t)u(r , t) = 0). The mean-field potential and related tensors can be 
decomposed into nuclear and Coulomb terms V(x,t) = V nU ci(x,t) + V CO ui(x 7 t). A nuclear potential can be derived 
as the first variation of the short-range interaction functional (as usual effective density dependent Skyrme forces) on 
density V nuc i(x, t) = 6£[n]/5n, This gives a nuclear potential as the function of the density n(r,t) = nor)(r,t). The 
incompressibility of the nuclear matter implies no = (2p^)/ (3nh 3 ) where no, pf being respectively, the density of 
nuclear matter and the Fermi momentum. Or, keeping in mind that an effective mean-field potential approximately 
constant inside the nuclei and has a sharp coordinate dependence in the surface region, to replace it by the surface 
term jll]] 'ygdivn, where fi — Vr/|Vr| is the unit outward normal on S, and the surface tension 7 S is related to mass 
formula coefficient via b s — Anr^ s ~ 22MeV. In the both cases the contour H(z,i) defines completely a nuclear 
potential on a nuclear surface. 
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V coul (x,t)=(jey m J dx' 

2 f 2 m „ /-SO',*) 



\x — X' 



tiq / dz' / do'o' I(cr,cr' , z'), 



I(a,a',z')= f 
Jo 





2 - ^ 



# (7) 

y/<7 2 + a' 2 — 2(7 cr' COS 



Concerning pressure tensors, we will not solve the kinetic equation (|l|), but will use that, this simplest picture can be 
reproduced by assuming the following factorised form of the distribution functions f(r,p,t) 

f°(r,p)^(27Th)- 3 V (r) ( f ) (p) f(r,p,t)^(27Th)- 3 r ] (r,t)^r,p,t), (8) 

where r)(r,t) = 9(T(r,t)) {r) a (r) = 8(T°(r))) is the step-function which is equal to one inside the surface T(r,t) 
(r°(r))and zero outside it. The distribution function f°(r,p) describes the equilibrium spherical symmetrical state 
(r°(r)) = a — \/i?o z 2 ), where Ro is the radius of a sphere of an equivalent volume (4i?Q = 3/ dzY,(z, t) 2 ). 
The distribution function f(r,p,t) corresponds to the dynamical picture. Both distribution function describe a 
homogenous distribution of a nuclear matter within the volume bounded by a sharp surface. The momentum- 
dependent part of the distribution functions could be parametrized as follows : 

4>°(p) = 9(j& -P 2 ), ^r,p,t) = 9(p 2 F -J2(^ + oiij(t))HQj)- (9) 

The forms chosen for the momentum-dependent part of f(r,p, t), (f°(r,p, t)), insures that it yields the current density 
with the collective velocity u(r, t)(u° =0). Of course, symmetry of the problem is defined by the symmetry of the all 
effective forces and an initial state. In a fluid at rest, only normal stresses are exerted, the normal stress is independent 
of the direction of the normal to the surface element across which it acts, and the equilibrium pressure tensor has 
the spherically symmetrical form via symmetry of the equilibrium distribution function P^(r, t) = (2 / 5)eFV ( r ) n o^ij ■ 
There no reason to expect these results to be valid for a fluid in motion, the tangential stresses are non-zero, in 
general, and the normal component of the stress acting across a surface element depends on the direction of the 
normal to the element. The simple notion of a pressure acting equally in all directions is lost in most cases of a fluid 
in motion. A deformations in the cartesian space define deformation in the momentum space and vice versa. We 
consider the only motion with 11 = and will assume the diagonal tensor aij. The dynamical pressure tensors are 
Pij(r, t) — (2/5)(l + au)~ 1 eFnor](r, t)6ij These equations show that the introduction of the functions etu{t) into the 
expression for <f>(r,p,t) accounts for the deformation of the Fermi-surface. 

The main hydrodynamical set of equations can be recast in the following form in the cylindrical coordinates 

du a du a du a 

-5T + U <r-H Vu z —- = F a , 

ot do oz 

du z du z du z 

-TTT + U °-B ^ u z-^~ = F z> 

ot OCT oz 

^ + ^ + ^ = 0, (10) 
da oz a 

where u a , u z {F„, F z ) are the projections of the velocity (forces) along the axis a and z respectively. 
The boundary conditions are 

hZh = - 7s (i- + — ) , fiEt = 0. (11) 

Here S is the stress tensor, which can be calculated using all the formulas above for the pressure tensor and the 
mean-field potentials. C\ and Ci are the principal radii of curvature. The direct evaluation gives 

ch + cl ~ e(i + £2)i/2 ~JT+ £2)3/2 ' *=a7' "-^Ih 2 (1 > 

e a -Y, z e z i Z z e a + e z . n + £ 2 t A -£ 2 n + t 

n = — i t — — ■ e a = — ■ e z — — , (13 

v/TTsf' VTTsf v/TTsf VT+^ Z 
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The final integro-differential equation of motion for a contour accumulates the coupled set of equations (|10[jl 

— jy t — = — -Q t N 2 (S,z,t) — — — = u a (S,z,t) (14) 

The evolution of the surface T(a, z, t) as a motion of a 2D space curve in the (a, z)-plane is a particular choice of a 
general 3D curve dynamics 

u = u t t + u n h + Ubb (15) 

where the tangent i, normal n, and binormal b = i x fi comprise the Frenet-Serret triad [Q, obeying the equations 
of differential geometry 

~ _ dr dt d 2 r _ „ _ d 2 r 

ds ' ds ds 2 ds 2 
db „ ofdr d 2 r\ d 3 r dn . „ s 

dS = - Tn > T = « 5 + ( 16 ) 

where k is the curvature, r is the torsion of the curve at arclength position s and time i. The direct evaluation gives 

^ = v /lTs!, r = 0, k= " 
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Axial symmetry of the problem reduces evolution to (n, t)-plane and (z <-> s) and (S(z, t) <-> k(s, t)) equivalence. This 
describes evolution of an axisymmctric surface pure geometrically in terms of the curvature and the arclength. Integro- 
differential form of Eqs. ( |10| - |l4| ) leads to strong nonlocality as a consequence of long-range part of the interaction. 



III. SUMMARY 



We consider an uniformly charged incompressible nuclear fluid bounded by a closed surface. It is shown that an 
evolution of an axisymmetric surface T(r, t) = a — S(z, t) = 0, r = (er, <f>, z) can be approximately reduced to a 
motion of a curve in the (er, z)-plane. A nonlinear integro-differential equation for the contour E(z,i) is derived. It 
is shown the direct correspondence between E(z, t) and a local curvature, that gives possibility to use methods of 
differential geometry to analyze an evolution of an axisymmetric nuclear surface. 

The goal of this short report is to present the principal line only. The extended paper is in progress. 

Work supported in part by Russian Foundation for Basic Research and Deutsche Forschungsgemeinschaft. 
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